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ABSTRACT 

Rosgen  stream  classification  and  rapid 
bioassessment  data  were  used  to  identify 
relationships  between  stream  reach  morphology 
and  the  macroinvertebrate  communities  of  high- 
quality,  cold  water  fisheries  (HQ-CWFs)  in  the 
Pennsylvania  portion  of  the  Susquehanna  River 
Basin.  Channel  substrate  composition,  channel 
entrenchment,  and  surface  water  slope  were  the 
stream  reach  morphology  variables  with  the 
greatest  influence  on  bioassessment  metric  scores. 
The  data  also  were  used  to  evaluate  the 
effectiveness  of  stream  order  (Strahler  method) 
and  Rosgen  stream  type  designation  as  stream 
reach  classification  variables.  Rosgen  stream 
classification  designation  was  more  effective  than 
stream  order  at  partitioning  biological  variation 
among  HQ-CWFs.  Taxonomic  richness  and 
modified  percent  mayfly  index  scores  were 
significantly  different  between  Rosgen  type  B and 
C streams.  A significant  difference  also  was 
observed  between  the  modified  Hilsenhoff  index 
scores  of  B3  and  B4  streams.  No  significant 
relationships  existed  between  stream  order  and 
bioassessment  metric  scores.  Rosgen 

classification  information  also  was  useful  for 
evaluating  the  overall  physical  habitat  conditions 
of  HQ-CWFs  and  identifying  streams  that  had 
deviated  from  stable,  '‘natural”  morphological 
conditions. 


INTRODUCTION 

In  an  effort  to  comply  with  federal 
antidegradation  requirements,  the  Pennsylvania 
Department  of  Environmental  Protection  (Pa. 
DEP)  developed  the  Pennsylvania  Special 


Protection  Waters  Program  (Shertzer,  1992).  The 
program,  which  includes  streams  or  watersheds 
that  exhibit  exceptional  water  quality  and  other 
environmental  features,  protects  selected 
waterways  by  establishing  higher  water  quality 
standards  than  those  that  typically  apply  to  the 
surface  waters  of  the  commonwealth.  Streams 
and  rivers  that  are  being  considered  for  inclusion 
in  the  program  usually  undergo  an  assessment  of 
their  "ecological  significance",  which  is 
determined  based  on  their  EPA  Rapid 
Bioassessment  Protocol  (RBP)  III  Score,  relative 
to  the  score  of  a comparable  reference  site. 
Reference  sites  represent  benchmarks  of  the  best 
attainable  conditions  that  can  be  expected  for 
similar  streams  within  a particular  ecological 
region  (Barbour  and  others,  1996).  Thus,  the 
delineation  of  ecological  regions  and  the 
identification  of  similar  streams  and  appropriate 
reference  sites  are  critical  components  of 
bioassessment  methodology. 

Barbour  and  others  (1996)  describe  the 
process  of  characterizing  regional  reference 
conditions  as  consisting  of  three  steps; 

1.  Classification  of  streams  based  on 
geographic  region  and  stream  fype; 

2.  Selection  of  the  best  available  sites  within 
a given  class  as  candidate  reference  sites; 
and 

3.  Characterization  of  reference  conditions 
based  on  biological  survey  data. 

These  authors  state  that  the  objectives  of  the  first 
step  in  this  process,  the  classification  of  streams 
based  on  geographic  region  and  stream  t\pe,  are 
to  reduce  the  complexity  of  biological  information 
and  improve  the  resolution  or  sensitivity  of 
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biological  surveys  by  partitioning  or  accounting 
for  variation  among  sites.  Furthermore,  they 
suggest  that  “the  best  classification  variables  are 
those  that  are  readily  obtained  from  maps  or 
regional  water  characteristics  such  as  ecoregion, 
gradient,  alkalinity,  and  hardness.”  Stream 
characteristics  that  are  readily  affected  by  human 
activities  or  occur  as  a biological  response  to 
physical  or  chemical  conditions  (i.e.,  land  use, 
habitat  condition,  or  nutrient  concentrations) 
should  not  be  used  as  classification  variables. 

Barbour’s  methodology  indicates  that  stream 
classification  should  include  variables  that 
describe  channel  features  at  a variety  of  different 
spatial  scales  such  as  geographic  region  and 
stream  type.  Gregory  and  others  (1991)  described 
the  spatial  hierarchy  of  drainage  basins  as 
consisting  of  channel  features  from  five  spatial 
scales:  (1)  network;  (2)  segment;  (3)  reach; 

(4)  channel  unit;  and  (5)  channel  subunit 
(Table  1).  However,  the  actual  number  and 
definition  of  spatial  scales  utilized  in  stream 
classification  programs  vary  and  may  be  dictated 
by  the  specific  objectives  of  the  program  and 


budgetary  and  manpower  constraints.  For 
example,  the  U.S.  Geological  Survey’s  National 
Water-Quality  Assessment  (NAWQA)  Program 
utilizes  a hierarchical  stream  habitat  assessment 
design  that  includes  characterizing  habitat  at  four 
spatial  scales:  (i)  basin;  (2)  stream  segment; 
(3)  stream  reach;  and  (4)  microhabitat. 

Currently,  the  Pa.  DEP  classifies  streams 
based  on  ecoregion  (or  subregion)  and  stream 
order  when  conducting  bioassessments  and 
ecological  significance  surveys  of  special 
protection  candidate  streams  (Shertzer,  1992). 
Biological  (macroinvertebrate)  data  collected 
during  these  surveys  is  limited  to  riffle  areas.  The 
classification  and  sampling  techniques  used  by  the 
Pa.  DEP  effectively  incorporate  classification 
variables  at  the  segment  and  channel  unit  spatial 
scales,  as  described  in  Table  1,  by  including 
ecoregion  designation  and  by  limiting  biological 
data  collection  to  riffles.  However,  stream  order, 
which  is  used  to  characterize  streams  at  the  reach 
spatial  scale,  does  not  characterize  stream  reach 
morphology,  which  is  a critical  component  of 
stream  ecosystems.  Stream  reach  morphology 


Table  1.  Spatial  Scales  Included  in  the  Hierarchical  Organization  of  Drainage  Basins  (Described  by 
Gregory  and  others,  1991) 


Spatial  Scale 

Number  of 
Channel  Widths 

Delineation  Features 

Network 

100,000 

Network  of  fluvial  and  valley  floor  landforms  within  a 
basin  extending  from  headwaters  to  large  rivers  or 
estuaries 

Segment 

1,000  - 10,000 

Regional  landform  patterns  that  reflect  different  agents  of 
landscape  formation  (scale  of  10  to  100  Km),  major 
topographic  discontinuities  such  as  high  gradient 
mountain  rivers  or  low  gradient  lowland  rivers 

Reach 

10  - 100 

Type  and  degree  of  local  constraint  imposed  on  the 
chWel  and  valley  floor  by  geomorphic  features 

Channel  Unit 

1 

Water-surface  slope,  channel  width/depth  ratio,  extent  of 
turbulent,  high-velocity  flow;  typically  classified  as 
riffles,  pools,  and  runs 

Channel  Subunits 

0.1 

Local  hydrauhc  features,  at  scales  less  than  the  active 
channel  width,  that  are  created  by  boulders,  logs,  or  gravel 
bars  and  exist  at  flows  less  than  those  required  to  fill  the 
active  charmel 

2 


influences:  (1)  Instream  habitat  conditions  and 
hydrograph  characteristics;  (2)  the  fluvial 
development  of  valley  floor  landforms  such  as 
flood  plains,  natural  levees,  and  terraces;  (3) 
riparian  vegetative  communities  and  the  quantity 
and  quality  of  organic  matter  entering  streams; 
and  ultimately;  (4)  aquatic  biological  community 
structure.  For  reviews  of  the  influence  of 
morphological  conditions  on  stream  and  riparian 
ecosystems,  see  Vannote  and  others  (1980)  and 
Gregory  and  others  (1991). 

The  Rosgen  Stream  Classification  System 
(Rosgen,  1996),  which  is  commonly  used  by  U.S. 
Department  of  Agriculture  (USD A),  Forest 
Service  hydrologists,  provides  a means  by  which 
morphologically  similar  stream  reaches  can  be 
identified.  This  system  utilizes  morphological 
variables  that  allow  for  the  classification  and 
assessment  of  stream  reach  morphology  at  four 
different  levels  of  resolution.  Rosgen  level  I and 
II  inventories  are  designed  for  classifying  stream 
reaches.  Level  I classification  is  a coarse-scale 


geomorphic  characterization  that  can  be 
conducted  using  information  taken  from 
topographic  and  landform  maps  and  from  aenal 
photographs.  Level  II  classification,  based  on 
field  data,  is  a more  detailed  morphological 
description  that  provides  a quantitative 
classification  of  stream  reaches.  Rosgen  level  I 
and  II  inventories  are  summarized  in  Table  2. 

Rosgen  level  III  and  lY  inventories  are  veiy 
detailed  and  build  upon  the  information  gathered 
during  level  I-II  inventories.  Level  III  inventories 
are  designed  for  assessing  the  existing  condition 
or  the  “state”  of  specific  stream  reaches.  Level  IV 
inventories  are  used  for  establishing  empirical 
relationships  pertaining  to  the  prediction  of 
Manning’s  “n”  values  and  studying  sediment 
transport  processes  and  hydraulic  geometiy- 
characteristics  of  streams  based  on  stream  t\pe  (as 
determined  from  level  11  data)  and  stream  state 
(based  on  level  III  data).  These  relationships  can 
then  be  extrapolated  to  similar  stream  reaches  for 
which  level  IV  data  are  not  available. 


Table  2.  Summary  of  Rosgen  Level  I and  II  Inventories 


Rosgen  Level  I (Geomorphic  Characterization) 

Primary  Information  Sources: 

Topographic  and  landform  maps 

Aerial  photography 

Delineation  Criteria; 

Channel  slope  (valley  slope/sinuosity) 

Channel  shape  (narrow/deep  vs.  wide/shallow) 

Channel  patterns  (single  thread,  multiple  thread,  anastamosed) 

Stream  Reach  Classes; 

Stream  types  "A"  through  "G" 

Rosgen  Level  II  (Morphological  Description) 

Primaiy'  Information  Sources: 

Field  Data 

Delineation  Criteria; 

Entrenchment  ratio 

Width/depth  ratio 

Sinuosity 

Channel  slope 

Chaimel  materials 

Stream  Reach  Classes: 

Stream  types  "Ala...A6a"  through  "Glc...G6c" 

To  accurately  classify  candidate  and  reference 
streams  for  the  purpose  of  comparing  their 
macro  invertebrate  communities,  it  is  essential  that 
these  streams  be  similar  at  a variety  of  spatial 
scales.  The  author  hypothesizes  that  the  Pa.  DEP 
could  strengthen  the  procedure  that  is  currently 
used  for  identifying  comparable  Special 
Protection  candidate  and  reference  sites  by 
incorporating  Rosgen  level  I or  II  stream  reach 
classification  designation  into  this  procedure.  The 
primary  objectives  of  this  project  are  to  evaluate 
the  potential  utility  of  the  Rosgen  classification 
system  as  a method  for  characterizing  HQ-CWFs 
at  the  stream  reach  spatial  scale  and  to  identify 
relationships  between  benthic  macroinvertebrate 
community  structure  and  stream  reach 
morphology  in  HQ-CWFs. 

METHODS 

Field  and  Laboratory  Methods 

A total  of  35  sample  sites  located  on  HQ- 
CWFs  were  sampled  between  September  1 1 and 
November  11,  1996  (Figure  1 and  Table  3). 
Sample  sites  were  located  in  the  Pennsylvania 
portion  of  the  Susquehanna  River  Basin  and 
represented  six  ecoregions  and  15  subregions 
(Figure  2).  At  each  sample  site,  stream  reach 
morphology  and  chemical  water  quality 
conditions  were  documented  in  the  field,  and 
benthic  macro  invertebrate  samples  were  collected 
for  analysis  in  the  laboratory. 

Chemical  water  quality 

Field  water  quality  measurements  included 
water  temperature,  dissolved  oxygen, 
conductivity,  pH,  alkalinity,  and  acidity. 
Dissolved  oxygen  was  measured  using  a YSI 
Model  55  dissolved  oxygen  meter.  Conductivity 
and  pH  were  measured  using  Cole-Parmer  Model 
1481-61  and  Model  5996-70  meters,  respectively. 
Alkalinity  was  measured  by  titrating  a known 
volume  of  sample  water  to  pH  4.5  with  0.02N 
H2SO4.  Acidity  was  measured  by  titrating  a 
known  volume  of  sample  water  to  pH  8.3  with 
0.02NNaOH. 


Benthic  macroinvertebrate  sampling 

Macroinvertebrate  communities  were  sampled 
using  a 0.3  m D-ffame  kick  net  to  collect 
organisms  dislodged  from  riffle  areas  by  physical 
agitation  of  approximately  0.3  m^  of  the 
streambed.  Two  samples  were  collected  at  each 
site,  one  from  an  area  of  high  velocity  and  the 
other  from  an  area  of  lower  velocity.  The  two 
samples  were  composited  and  preserved  in  a 
solution  of  isopropyl  alcohol  and  glycerin  for 
laboratory  analysis.  In  the  laboratory, 
100-organism  subsamples  were  taken  from 
composite  samples  using  a gridded  pan  and  a 
random  numbers  table.  The  organisms  contained 
in  the  subsamples  were  identified  to  genus  (except 
for  Chironomidae  and  Simuhidae)  using 
taxonomic  keys  by  Merrit  and  Cummins  (1996) 
and  Peckarsky  and  others  (1990)  and  enumerated. 

Stream  reach  morphology 

At  each  sample  site  (sample  reach),  stream 
reach  morphological  conditions  were  character- 
ized using  Rosgen  level  II  classification  data 
(Rosgen,  1996).  Level  n classification  data 
consisted  of  field  measurements  of  channel  cross- 
section  dimensions,  sinuosity,  slope,  and  the  size 
of  channel  materials.  Channel  cross-section 
dimensions,  sinuosity,  and  slope  were  determined 
using  a surveying  rod,  level,  and  tape  measure.  At 
each  site,  sinuosity,  slope,  and  channel  materials 
were  assessed  over  the  entire  sample  reach. 
Sample  reaches  were  approximately  equal  to  20 
bankfull  channel  widths  in  length.  Sample  reach 
end-points  were  delineated  such  that  reaches 
started  and  ended  in  the  same  type  of  channel  unit 
(e.g.,  head-of-riffle  to  head-of-riffle). 

At  each  sample  site,  a stream  channel  cross- 
section,  located  in  a relatively  straight, 
representative  segment  of  the  sample  reach,  was 
surveyed.  While  surveying  cross-sections, 
bankfull  stage  was  identified  using  field  stage 
indicators  such  as  those  described  by  Dunne  and 
Leopold  (1978),  Harrelson  and  others  (1994),  and 
Rosgen  (1996).  Some  of  these  indicators 
included:  floodplain  elevations,  the  tops  of 
depositional  features  such  as  point  bars,  breaks  in 
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Figure  1.  Map  Showing  the  Location  of  Sample  Sites  in  the  Susquehanna  River  Basin 


Table  3.  Descriptions  of  35  Sample  Site  Locations 
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Figure  2.  Map  Showing  the  Ecoregions  of  Pennsylvania  and  the  Location  of  Sample  Sites 


streambank  slope,  changes  in  streambank  particle 
size,  and  staining  of  rocks.  Cross-section  data 
were  analyzed  using  XSPRO  stream  charmel 
cross-section  analysis  software  (Grant  and  others, 
1992)  to  obtain  cross-section  dimensions. 

Entrenchment  and  width/depth  ratios  were 
derived  from  channel  cross-section  data. 
Entrenchment  ratio  describes  the  vertical 
containment  of  the  stream  channel  (the  degree  to 
which  it  is  incised  in  the  valley  floor). 
Entrenchment  ratio  was  calculated  as  the  ratio  of 
the  width  of  the  flood-prone  area  to  the  surface 
width  of  the  bankfull  channel.  The  flood-prone 
area  width  was  measured  at  the  elevation  that 
corresponded  to  twice  the  maximum  depth  of  the 
bankfull  charmel  (Figure  3).  Width/depth  ratio 
was  calculated  as  the  ratio  of  the  bankfull  surface 
width  to  the  mean  depth  of  the  bankfull  charmel. 

Sinuosity  was  calculated  as  the  ratio  of  stream 
charmel  length  to  valley  length  of  the  study  reach. 
Stream  length  was  determined  in  the  field  by 
measuring  the  length  of  each  charmel  unit  (riffle, 
run,  or  pool)  in  the  study  reach  using  a tape 
measure  extended  in  alignment  with  the  charmel 
thalweg.  Valley  length  was  determined  by 
measuring  the  straight  hne  distance  between  the 
upstream  and  downstream  end-points  of  the  reach, 
which  were  surveyed  in  the  field  and  plotted  in 
the  office.  Water  surface  slope  was  calculated  as 
the  change  in  surface  water  elevation  per  unit 
stream  length.  Surface  water  elevation  at  the 
upstream  and  downstream  end-points  of  each 
study  reach  were  surveyed  in  the  field.  Stream 
length  was  determined  in  the  field  using  a tape 
measure,  as  described  above. 

Dominant  charmel  materials  were  determined 
based  on  the  particle  size  distribution  of  bankfull 
channel  materials  using  the  modified  Wolman 
“pebble-count”  method,  as  described  by  Rosgen 
(1996).  The  intermediate  axis  diameter  of 
approximately  100  randomly  selected  substrate 
particles  was  measured  and  tabulated  in  the  field. 
The  number  of  particles  selected  from  a given 
channel  unit  (riffle,  run,  or  pool)  was  based  on  the 
charmel  unit  configuration  of  the  stream  reach, 
such  that  charmel  units  were  sampled  in 


proportion  to  their  occurrence  in  the  reach.  The 
dominant  particle  size  of  the  stream  reach  was 
determined  as  the  indexed  dso  value  from  a 
cumulative  percent  plot  of  the  pebble  count  data 
(Appendix  A).  Rosgen  level  11  stream  type 
designations  were  determined  using  the  criteria 
summarized  in  Table  4. 

Ecoreqion  designation,  stream  order,  and 

watershed  characteristics 

Sample  site  ecoregion  (subregion)  designation 
and  watershed  characteristics  (e.g.,  drainage  area, 
percent  ecoregion  composition)  were  determined 
using  PC  ArcView  spatial  analyst  geographic 
information  system  (GIS)  software  and  readily 
available  geographic  data.  Stream  order  (Strahler 
method)  was  delineated  from  United  States 
Geological  Survey  7.5  minute  quadrangle  maps. 
Intermittent  streams  were  not  included  in  stream 
order  designations. 

Data  Analysis 

All  statistical  tests  and  simple  linear 
correlation  analyses  were  conducted  in 
accordance  with  Zar  (1984)  using  MINITAB 
statistical  software  (MINITAB,  1996).  Box-and- 
whisker  plots  were  generated  using  this  same 
software.  Significant  results  were  reported  at  the 
95  percent  confidence  level  (p<0.05).  When 
statistically  significant  analysis  of  variance 
(ANOVA)  results  were  obtained,  differences 
among  the  sample  groups  were  evaluated  using 
Tukey  multiple  comparison  tests  and  box-and 
whisker  plots.  Box-and- whisker  plots  graphically 
display  the  median,  minimum,  maximum,  25^^  and 
75*^  percentiles  (also  referred  to  as  the  T‘  and  3^'^ 
quartiles,  respectively),  and  outliers  of  a given 
sample  group.  Box-and- whisker  outliers  are 
points  outside  of  the  lower  and  upper  limits  of  the 
sample  group.  Sample  group  limits  were 
determined  as  follows: 

Lower  Limit:  T'  quartile  - 1.5  (3'^*  quartile  - T' 
quartile) 

Upper  Limit:  3'“^  quartile  - 1.5  (3'‘^  quartile  - U‘ 
quartile) 
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Figure  3.  Illustration  of  Entrenchment  Ratio  Calculated  From  Bankfull  and  Flood-Prone  Channel 
Widths 
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Table  4.  Key  Used  for  Rosgen  Level  I and  II  Classification  of  Sample  Sites 
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Values  of  Entrenchment  Ratio  and  Sinuosity  can  vary  by  +/-0.2  units;  values  for  Width/Depth  Ratio  can  vary  by+/-  2.0  units.  Source:  Modified  from  Rosgen  (1996). 


Detrended  correspondence  analysis  (DCA) 
was  used  to  reveal  outliers  in  the  data  set  based  on 
macroinvertebrate  taxa  abundance.  Gauch  (1982) 
describes  outliers  as  peculiar  samples  whose 
relationships  to  other  samples  in  the  data  set  are 
not  expressed  by  the  information  in  the  data. 
Furthermore,  this  author  states  that,  in  general, 
outliers  tend  to  compress  the  other  samples  in  the 
data  set  into  an  undesirably  small  area  of  the 
ordination  and  obscure  the  analysis  of  the  data  set 
as  a whole,  and  thus,  should  be  deleted  from  the 
data  set.  Since  the  focus  of  this  study  is  on  the 
relationships  between  stream  reach  morphology 
and  the  benthic  macroinvertebrate  communities  of 
HQ-CWFs,  samples  that  were  primarily 
influenced  by  water  quality  conditions  or  that  had 
obviously  degraded  channel  morphology  or 
macroinvertebrate  communities  were  removed 
from  the  analysis.  DCA  ordinations  were 
produced  using  software  developed  by  Kovach 
(1993).  Rare  taxa,  those  present  at  less  than 
5 percent  of  the  sites,  were  removed  prior  to 
miming  ordination  analyses  (Gauch,  1982). 

Five  metrics  were  used  to  describe  sample  site 
benthic  macroinvertebrate  communities.  These 
metrics,  which  included  taxa  richness,  modified 


EPT  index,  modified  Hilsenhoff  index,  percent 
dominant  taxa,  and  modified  percent  ma\flv 
index,  are  summarized  in  Table  5.  The  Pa.  DEP 
was  considering  this  suite  of  metrics  for  assessmg 
the  ecological  significance  of  Special  Protection 
candidate-streams  at  the  start  of  this  project 
(Pa.  DEP,  Unpublished). 

RESULTS 

Data  Summary  and  the  Identification  of 
Outliers 

Rosgen  level  I classification  delineates  the 
35  sample  sites  included  in  this  study  into  five 
different  stream  types.  The  majority  of  the 
sample  sites  are  either  B or  C stream  types,  which 
together  account  for  80  percent  of  the  sites 
included  in  the  study.  Rosgen  level  II  stream 
classification  further  delineates  the  35  sample 
sites  into  13  different  stream  types  consisting 
primarily  of  types  B4c  (20  percent),  C3 
(14  percent),  C4b  (11  percent),  and  C4 
(11  percent)  (Figure  4).  Rosgen  stream  reach 
classification  information  is  summarized  in 
Table  6.  Stream  channel  morphology  and  water 
quality  data  are  listed  in  Appendixes  A and  B. 


Table  5.  Summary  of  the  Bioassessment  Metrics  Used  to  Describe  Sample  Site  Benthic 

Macroinvertebrate  Communities 


Metric 

Description 

Taxa  Richness 

Total  number  of  different  taxa. 

Modified  EPT  Index 

Number  of  pollution  intolerant  (Hilsenhoff  tolerance  value  < 5) 
Epheraeroptera  (mayfly),  Plecoptera  (stonefly),  and  Trichoptera  (caddisflv) 
taxa 

Modified  Hilsenhoff  Index 

Hilsenhoff  index  of  pollution  tolerance  with  pollution  tolerance  values 
modified  for  Pennsylvania 

Percent  Dominant  Taxa 

Relative  abundance  of  the  dominant  taxa  expressed  as  percent  of  the  100- 
organism  subsample 

Modified  Percent  Ivlayfly  Index 

Relative  abundance  of  pollution  intolerant  (Hilsenhoff  tolerance  value  < 5) 
mayflies  e.xpressed  as  percent  of  the  100-organism  subsample 
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Figure  4.  Histograms  of  Rosgen  Level  I and  II  Classification  Results  of  35  Sample  Sites 
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Table  6.  Summary  of  Rosgen  Stream  Reach  Classification  Data  of  35  Sample  Sites 


Sample 

Site 

Entrenchment 

W/D 

Ratio 

Sinuosity 

Slope 

Channel 

Material 

(dso) 

Rosgen 
Level  1 

Class 

Rosgen 
Level  II 

Class 

Stream  j 

Order 

1:24,000 

ASPH 

1.8 

28 

1.1 

0.020 

Cobble 

B 

B3 

3 
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1.4 

30 

1.1 

0.029 

Cobble 

B 

B3 

2 

CORK 

2.0 

22 

1.2 

0.019 

Gravel 

B 

B4 

1 
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1.4 

25 

1.1 

0.036 
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B 

B3 

1 

DWMN 

4.5 

22 

1.4 

0.029 

Gravel 

C 

C4b 

2 
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1.7 

13 

1.0 

0.002 

Gravel 

B 

B4c 

3 

EBLR 

3.8 

16 

1.0 

0.010 

Cobble 

C 

C3 

2 

EBRB 

1.3 

17 

1.0 

0.019 

Cobble 

B 

B3c 

2 

GAZZ 

1.4 

26 

1.4 

0.008 

Gravel 

B 

B4c 

3 

GRLE 

4.5 

14 

1.4 

0.014 

Cobble 

C 

C3 

1 

GRNL 

1.3 

27 

1.5 

0.014 

Gravel 

B 

B4c 

2 

GRNU 

1.3 

31 

1.1 

0.012 

Gravel 

B 

B4c 

2 

HAML 

1.3 

20 

1.1 

0.003 

Gravel 

F 

F4 

3 

HAMU 

1.7 

11 

1.0 

0.003 

Sand 

G 

G5c 

2 

HARV 

2.0 

23 

1.2 

0.005 

Gravel 

B 

B4c 

2 

HORS 

4.5 

8 

1.2 

0.025 

Sand 

E 

E5b 

1 

HUTT 

2.8 

14 

1.1 

0.040 

Gravel 

C 

C4b 

2 

IRON 

4.5 

21 

1.1 

0.019 

Gravel 

C 

C4 

2 

LEBO 

2.9 

25 

1.2 

0.015 

Cobble 

c 

C3 

3 

LMOR 

2.0 

17 

1.2 

0.027 

Gravel 

B 

B4 

2 

LTRT 

1.7 

14 

1.1 

0.006 

Sand 

G 

G5c 

1 

MNTN 

4.5 

14 

1.1 

0.005 

Gravel 

C 

C4 

3 

MURY 

1.9 

16 

1.0 

0.014 

Gravel 

B 

B4c 

1 

PEQU 

1.6 

16 

1.0 

0.017 

Gravel 

B 

B4c 

2 

RICE 

2.3 

14 

1.2 

0.040 

Gravel 

C 

C4b 

1 

SIKM 

2.8 

21 

1.0 

0.016 

Cobble 

C 

C3 

3 

STSN 

4.5 

15 

1.2 

0.009 

Cobble 

C 

C3 

3 

TIPT 

1.6 

26 

1.1 

0.010 

Cobble 

B 

B3c 

2 

TOML 

4.5 

7 

1.0 

0.001 

Gravel 

E 

E4 

1 

TOMU 

2.4 

33 

1.1 

0.020 

Gravel 

C 

C4b 

1 

UEBC 

2.0 

8 

1.2 

0.014 

Gravel 

G 

G4c 

2 

UHAR 

1.4 

23 

1.1 

0.041 

Cobble 

B 

B3a 

1 

URRB 

3.6 

17 

1.1 

0.018 

Gravel 

C 

C4 

2 

WILL 

4.5 

18 

1.2 

0.005 

Gravel 

C 

C4 

1 

WTDR 

4.5 

11 

1.3 

0.004 

Gravel 

E 

E4 

2 
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A total  of  83  taxa  from  43  families  were 
included  in  the  macroinvertebrate  data  collected  at 
the  35  sample  sites  (Appendix  C).  Figure  5 shows 
the  DCA  ordination  of  sample  sites  based  on 
macroninvertebrate  taxa  abundance  (data  set  1). 
This  ordination  suggests  LTRT,  HAMU,  HAJVIL, 
WILL,  HORS,  HUTT,  and  HARV  may  be  outliers 
in  the  data  set.  Correlation  analysis  of  channel 
morphology  and  water  quality  data  with  sample 
site  DCA  ordination  scores  suggests  that 
macroinvertebrate  community  structure  is 
primarily  influenced  by  conductivity  (r  = 0.775, 
p < 0.001)  and  alkalinity  (r  = 0.737,  p < 0.001), 
and  to  a lesser  extent,  degraded  channel  substrate 
(Figure  6).  LTRT,  HAMU,  HAML,  WILL, 
HORS,  and  HUTT  were  determined  to  be  outliers 
based  on  the  rationale  described  in  Table  7. 
HARV  was  retained  since  it:  (1)  is  not  influenced 
by  limestone  geology;  (2)  is  characterized  by 
stable,  natural  stream  reach  morphology;  and 
(3)  supports  a healthy  macroinvertebrate 
community. 


A second  data  set  (data  set  2),  consisting  of 
29  sites,  was  constructed  by  removing  the  six 
outliers  mentioned  above  from  the  original  data 
set  of  35  sites.  DCA  performed  on  data  set  2 
produced  a DCA  ordination  of  29  sample  sites 
that  does  not  suggest  the  presence  of  any 
additional  outliers  (Figure  7).  Correlation 
analysis  of  DCA  axis  scores,  morphological 
characteristics,  and  water  quality  data  indicate  that 
macroinvertebrate  eommunity  structure  in  data  set 
2 is  predominantly  influenced  by  dissolved 
oxygen,  percent  cobble  substrate,  acidity,  and 
percent  riffle.  The  metrie  scores  of  sample  sites 
are  listed  in  Table  8.  The  metric  scores  of  data  set 
1 (35  sites)  and  data  set  2 (29  sites,  outliers 
identified  by  DCA  removed)  are  summarized  in 
Table  9 and  Figure  8. 
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Figure  5.  Detrended  Correspondence  Analysis  (DCA)  Ordination  of  35  Sample  Sites  Based  on 
Macroinvertebrate  Taxonomic  Abundance 
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Figure  6.  Scatter  Diagrams  of  DCA  Axr/s  1 Scores  Vs.  the  Water  Quality  and  Channel  Morphology  Parameters  Most  Strongly  Correlated 
With  the  DCA  Axis  1 Scores  of  35  Sample  Sites 


Table  7.  Summary  of  the  Six  Outlier  Sample  Sites  and  the  Justification  for  Their  Removal  From  the 
Original  Data  Set  of  35  Sample  Sites 


Site 

Justification  for  Removal  From  Data  Set 

LTRT 

Elevated  conductivity  and  alkalinity,  geochemically  influenced 

Site  is  located  in,  and  entire  watershed  consists  of,  Ecoregion  67a  (limestone  or  dolomite  valleys 
subregion) 

Degraded  charmel  morphology  (G5c  stream  type  in  disequilibrium) 

HANOI 

Elevated  conductivity  and  alkalinity,  geochemically  influenced 

Site  is  located  in,  and  51  percent  of  its  watershed  consists  of,  Ecoregion  67a  (hmestone  or  dolomite 
valleys  subregion) 

Degraded  channel  morphology  (G5c  stream  type  in  disequilibrium) 

HANIL 

Elevated  conductivity  and  alkalinity,  geochemically  influenced 

60  percent  of  the  site's  watershed  consists  of  Ecoregion  67a  (hmestone  or  dolomite  valleys 
subregion) 

Degraded  charmel  morphology  (F4  stream  type  in  disequiUbrium) 

WILL 

Elevated  conductivity  and  alkalinity,  geochemically  influenced 

12  percent  of  the  site's  watershed  consists  of  Ecoregion  67a  (hmestone  or  dolomite  valleys 
subregion) 

HORS 

Degraded  charmel  substrate  composition  (52  percent  sand) 

Poor  metric  scores 

HUTT 

Poor  metric  scores 
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Figure  7.  DC  A Ordination  of  Data  Set  2 Sample  Sites  (Data  set  2 consists  of  29  sites  and  was 

constructed  by  removing  the  six  outliers,  identified  by  DCA.  from  the  original  data  set  of 
35  sample  sites  (data  set  1).) 
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Table  8.  Metric  Scores  of  35  Sample  Sites 


Site 

Tax  Rich  Mod  EPT 

Mod  Hits 

% Dom 

Mod  % May 

ASPH 

25 

15 

2.37 

21 

21 

BALD 

19 

10 

2.28 

16 

15 

CORK 

16 

11 

1.66 

27 

36 

DEEP 

24 

15 

1.52 

25 

25 

DWMN 

7 

4 

2.81 

29 

0 

EBCC 

17 

10 

3.43 

29 

25 

EBLR 

20 

12 

2.48 

19 

44 

EBRB 

28 

16 

1.98 

13 

20 

GP3LZ 

25 

9 

3.75 

14 

10 

GRLE 

13 

5 

3.10 

24 

8 

GRNL 

22 

12 

3.00 

18 

15 

GRNU 

20 

10 

2.80 

16 

22 

HAML* 

10 

4 

4.28 

41 

10 

HAMU* 

10 

2 

5.52 

31 

1 

HARV 

16 

10 

3.32 

20 

42 

HORS* 

3 

0 

7.07 

79 

0 

HUTT* 

5 

3 

1.68 

48 

48 

IRON 

21 

10 

2.69 

15 

3 

LEBO 

19 

13 

2.91 

27 

12 

LMOR 

23 

12 

3.56 

16 

18 

LTRT* 

7 

0 

6.54 

56 

0 

MNTN 

21 

12 

3.58 

20 

16 

MURY 

20 

9 

3.00 

21 

17 

PEQU 

19 

8 

3.72 

49 

7 

RICE 

14 

11 

2.07 

22 

15 

SIXM 

17 

9 

2.70 

21 

3 

STSN 

17 

9 

2.38 

29 

2 

TIPT 

21 

11 

2.31 

37 

10 

TOML 

16 

7 

3.41 

12 

3 

TOMU 

17 

10 

2.77 

13 

2 

UEBC 

17 

7 

2.86 

22 

9 

UHAR 

15 

10 

3.08 

21 

23 

URRB 

20 

11 

2.41 

20 

14 

WILL* 

20 

10 

3.23 

27 

45 

WTDR 

22 

13 

2.44 

15 

5 

* Outlier 


Table  9.  Summary  of  Metric  Scores  of  Data  Sets  1 (35  Sample  Sites)  and  2 (29  Sample  Sites, 

6 Outliers  Identified  by  DCA  Removed) 


Taxonomic 

Richness 

Modified 

EPT 

Index 

Modified 

Hilsenhoff 

Index 

Per( 

Dom 

Ind 

:ent 

nant 

ex 

Modified  Percent 
Mayflly 

Index 

Data  Set 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

Number  of  Sites 

35 

29 

35 

29 

35 

29 

35 

29 

35 

29 

Minimum 

3 

7 

0 

4 

1.52 

1.52 

12 

12 

0 

0 

Maximum 

28 

28 

16 

16 

7.07 

3.75 

79 

49 

48 

44 

Mean 

17.3 

19.0 

9.1 

10.4 

3.10 

2.80 

26.1 

21.8 

15.6 

15.2 

Median 

19.0 

19.0 

10.0 

10.0 

2.90 

2.80 

21.0 

21.0 

14.0 

15.0 
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Figure  8.  Summary  of  Metric  Scores  of  Data  Sets  1 (35  Sample  Sites)  and  2 (29  Sample  Sites, 
6 Outliers  Identified  by  DCA  Removed) 
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Relationships  Between  Rosgen 
Classification  and  Metric  Scores  and 
Between  Stream  Order  and  Metric  Scores 

Twenty-six  of  the  29  sample  sites  included  in 
data  set  2 are  classified  as  Rosgen  level  I type  B 
or  C streams.  These  26  sites  were  used  to 
evaluate  the  effectiveness  of  stream  order  and 
Rosgen  stream  type  designation  as  stream  reach 
classification  variables  (data  set  3).  The  results  of 
statistical  analyses  to  test  for  significant 
differences  in  the  metric  scores  of  T‘  and  2"^^  order 
streams  vs.  3^^^  order  streams  and  between  type  B 
vs.  type  C streams  are  summarized  in  Table  10. 
These  analyses  indicate  that  there  is  a significant 
difference  in  the  mean  taxonomic  richness  and 
modified  percent  mayfly  metric  scores  of  Rosgen 
level  I type  B and  type  C streams  and  that  there  is 
no  significant  difference  in  any  of  the  metric 
scores  of  and  2"‘*  order  streams  vs.  S"'*  order 
streams.  Metric  scores  of  data  set  3 sample  sites, 
based  on  stream  order  and  level  I stream  type 
designation,  are  summarized  in  Figure  9. 

Rosgen  level  II  classification  refines  the 
major  stream  types  (types  A through  G)  into 
subcategories  based  on  channel  materials  and 
water  surface  slope.  For  example,  a level  11  B3c 


stream  reach  is  a B stream  type  with  cobble 
substrate  and  a surface  water  slope  of  less  than 
0.020.  Rosgen  level  II  classification  of  data  set  3 
sample  sites  separates  these  26  type  B and  C 
streams  into  eight  different  level  II  stream  types. 
Due  to  the  excessive  fragmentation  of  the  data  set 
and  the  low  number  of  streams  in  many  of  the 
level  II  stream  type  categories,  level  II 
classification  results  did  not  support  further 
investigation  into  the  relationships  between  level 
n stream  types  and  macroinvertebrate  community 
structure.  However,  excluding  the  surface  water 
slope  subcategory  designation  and  classifying 
sample  sites  based  on  major  stream  type  category 
and  channel  materials  partitions  the  26  type  B and 
C sites  into  one  of  four  “modified  level  II”  stream 
type  categories  (Figure  10). 

ANOVA  results  indicate  that  modified 
Hilsenhoff  index  is  the  only  metric  with  scores 
that  are  significantly  different  among  the  four 
modified  level  11  stream  types  (Table  11  and 
Figure  11).  Tukey  multiple  comparison  test 
results  indicate  that  mean  Hilsenhoff  scores  of  B3 
steams  are  significantly  different  than  those  of  B4 
streams,  but  that  there  is  no  significant  difference 
among  mean  Hilsenhoff  scores  of  B3,  C3,  and  C4 
or  among  B4,  C3,  and  C4  streams. 


TablelO.  Summary  of  F-Test  of  Sample  Variances  and  t-Tests  of  Sample  Means  to  Test  for 

Significant  Differences  in  the  Metric  Scores  of  Data  Set  3 Sample  Sites  Based  on  Stream 
Order  and  Rosgen  Level  I Stream  Type  Designation 


F-Test:  Two-Sample  for  Variances 

t-Test:  Two-Sample  Assuming  Equal  Variances 

Parameter 

Fcalculated 

Fcrttical 

tcalculated 

tcritlcal 

P(T<=t)  two-tail 

Type  B Type  C Streams 

Taxa  Richness 

1.22 

Fo.05(2),  10.14  =3.15 

2.366* 

to.05  (2),  24=  2.064 

0.026 

Modified  EPT 

1.40 

Fo.05(2),  10,14  = 3.15 

1.521 

to.05  (2),  24  = 2.064 

Modified  HUsenhoff  Index 

3.28 

Fo.05(2),  14,10  “ 3.55 

0.275 

to,05  (2),  24  = 2.064 

% Dominant  Taxa 

3.39 

Fo.05(2),  14,10  ~ 3.55 

0.354 

to  05  (2).  24  = 2.064 

Modified  % Mayfly 

1.77 

Fo.05(2),  10,14  =3.15 

2.237* 

to.05  (2).  24  = 2.064 

0.035 

1st  and  2nd  Order  ra:  3rd  Order  Streams 

Taxa  Richness 

1.62 

Fo.05(2),  18.6  = 5.20 

-0.751 

lo.05  (2),  24  = 2.064 

Modified  EPT 

1.38 

Fo.05(2),  18,6  ~ 5.20 

-0.530 

to.05(2),  24  = 2.064 

Modified  Hilsenhoff  Index 

1.10 

Fo.05(2).  18.6  ~ 5.20 

-1.357 

to.05  (2),  24  = 2.064 

% Dominant  Taxa 

2.50 

Fo.05(2),  18,6  5.20 

-0.235 

lo.05  (2),  24  = 2.064 

Modified  % Mayfly 

2.09 

Fo.05f2),  18,6  “ 5.20 

0.966 

lo.05('2),24  = 2.064 

* p<0.05 
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Figure  9.  Summary  of  Data  Set  3 Metric  Scores  Based  on  Stream  Order  and  Rosgen  Level  I Stream 
Type  Designation 
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Figure  9.  Summary  of  Data  Set  3 Metric  Scored  Based  on  Stream  Order  and  Rosgen  Level  I Stream 

Type  Designation — Continued 
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Table  11.  ANOVA  Results  of  Metric  Scores  of  Modified  Level  II  Stream  Types 


Source  of  Variation 

SS 

df 

MS 

F 

P-value 

F crit 

Single  Factor  ANOVA  of  Taxa  Richness  Scores 

Between  Groups 

108.157 

3 

36.052 

2.169 

0.120 

3.049 

Within  Groups 

365.689 

22 

16.622 

Total 

473.846 

25 

Single  Factor  ANOVA  of  Modified  EPT  Index  Scores 

Between  Groups 

42.206 

3 

14.069 

2.305 

0.105 

3.049 

Within  Groups 

134.256 

22 

6.103 

Total 

176.462 

25 

Single  Factor  ANOVA  of  Modified  Hilsenhoff  Index  Scores 

Between  Groups 

2.824 

3 

0.941 

3.295* 

0.039 

3.049 

Within  Groups 

6.283 

22 

0.286 

Total 

9.107 

25 

Single  Factor  ANOVA  of  Percent  Dominant  Taxa  Index  Scores 

Between  Groups 

60.487 

3 

20.162 

0.291 

0.831 

3.049 

Within  Groups 

1,523.667 

22 

69.258 

Total 

1,584.154 

25 

Single  Factor  ANOVA  of  Modified  Percent  Mayfliy  Index  Scores 

Between  Groups 

683.751 

3 

227.917 

1.863 

0.165 

3.049 

Within  Groups 

2,692.133 

22 

122.370 

Total 

3,375.885 

25 

* Significant  at  p<0.05 
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Modified  EPT  Index  vs.  Stream  Type 
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Figure  11.  Summary  of  Data  Set  3 Metric  Scores  Based  on  Modified  Level  II  Stream  Type 
Designation 
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Relationships  Between  Stream  Reach 
Channel  Morphology  and  Metric  Scores 

Channel  substrate  particle  size,  entrenchment, 
and  slope  are  the  stream  reach  morphology 
variables  that  have  the  greatest  influence  on  the 
metric  scores  of  the  26  HQ-CWFs  included  in 
data  set  3 (non-outlier,  type  B and  C streams). 
Taxa  richness,  modified  EPT,  and  modified 
Hilsenhoff  index  scores  are  significantly 
correlated  to  substrate  particle  size,  as  described 
by  di5,  d35,  and  d5o  values.  Modified  EPT  and 
modified  Hilsenhoff  index  scores  also  are 
significantly  related  to  the  amount  of  very  coarse 
(percent  boulder)  materials  present  in  the  stream 
channel.  In  addition,  modified  Hilsenhoff  index 
scores  are  influenced  by  amount  of  silt  and  clay 
(percent  silt/clay)  in  the  stream  channel.  Channel 
morphology  parameters  that  are  significantly 
correlated  with  the  metric  scores  of  data  set  3 
sample  sites  are  summarized  in  Table  12. 

In  addition  to  channel  substrate  composition, 
surface  water  slope  and  channel  entrenchment 
substantially  influence  the  metric  scores  of  sample 
sites.  A statistically  significant  correlation  exists 


between  surface  water  slope  and  modified 
Hilsenhoff  index  scores.  Entrenchment  ratio, 
which  is  the  primary  stream  reach  morphology 
variable  that  distinguishes  Rosgen  level  I type  B 
streams  from  type  C streams,  is  significantly 
correlated  to  taxa  richness  and  modified  EPT 
index  scores. 

Although  the  mean  Hilsenhoff  index  scores  of 
level  I type  B and  C streams  are  not  significantly 
different,  the  mean  Hilsenhoff  index  scores  of  the 
modified  level  II  stream  types  differ  significantly. 
Modified  level  II  data  indicate  that  stream  reaches 
characterized  by  low  surface  water  slope  values, 
in  conjunction  with  small  chaimel  substrate 
materials,  tend  to  have  relatively  high  modified 
Hilsenhoff  index  scores.  This  relationship  is 
particularly  evident  in  modified  level  n B4 
streams  (Figure  12).  ANOVA  and  Tukey  multiple 
comparison  tests  indicate  that  Hilsenhoff  index 
scores,  d35  particle  size,  and  surface  water  slope 
are  significantly  different  between  B3  and  B4 
stream  types,  and  the  d35  particle  size  of  B4  and 
C4  streams  are  different  than  those  of  B3  and  C3 
streams. 


Table  12.  Summary  of  Correlation  Analysis  of  Stream  Reach  Morphology  Parameters  and  Metric 
Scores  of  Data  Set  3 Sample  Sites  (T able  includes  morphology  parameters  that  are 
significantly  correlated  (p<  0. 05)  to  at  least  one  of  the  five  metrics.) 


Tax  Rich 

Mod  EPT 

Mod  Mils 

% Dom 

Mod  % May 

Entrenchment 

-0.429* 

-0.410* 

0.055 

0.026 

-0.326 

Slope 

-0.184 

0.089 

-0.413* 

-0.085 

1 

p 

o 

dl5 

0.437* 

0.579* 

-0.490* 

-0.142 

0.182 

d35 

0.376 

0.486* 

-0.518* 

0.040 

0.093 

d50 

0.400* 

0.531* 

-0.473* 

-0.074 

0.139 

Percent  Silt/Clay 

0.001 

-0.320 

0.428* 

-0.059 

0.083 

Percent  Boulder 

0.231 

0.398* 

-0.412* 

-0.076 

0.137 

* p<  0.05 
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Modified  Hilsenhoff  Index  vs.  Stream  Type 
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IVater  Slope,  and  Modified  Level  II  Stream  Type  Designation 


DISCUSSION 

Three  of  the  five  stream  reach  morphology 
parameters  used  in  the  Rosgen  classification 
system  are  significantly  correlated  with  at  least 
one  of  the  metrics  included  in  this  study.  As  d50 
particle  size  increases,  taxonomic  richness  and 
modified  EPT  index  scores  increase,  and  modified 
Hilsenhoff  index  scores  decrease.  As 
entrenchment  ratios  increase,  taxonomic  richness 
and  modified  EPT  index  scores  decrease.  Low 
surface  water  slope  values  correspond  to  elevated 
modified  Hilsenhoff  index  scores,  especially  in 
streams  where  low  surface  water  slope  values 
coincide  with  channel  substrates  that  are 
characterized  by  low  dss  particle  size. 

In  contrast  to  the  Rosgen  stream  reach 
classification  parameters,  stream  order 
designation  is  not  significantly  correlated  with  any 
of  the  five  metrics  included  in  this  study. 
Furthermore,  there  are  no  significant  differences 
between  the  mean  metric  scores  of  sample  sites 
when  they  are  grouped  based  on  stream  order 
designation.  However,  significant  differences 
exist  among  the  mean  scores  of  several  of  the 
metrics  when  sample  sites  are  grouped  based  on 
Rosgen  stream  type  designation.  Thus,  the  results 
of  this  study  indicate  Rosgen  stream  classification 
designation  is  more  effective  at  partitioning 
biological  variation  among  HQ-CWFs  than  stream 
order. 

Entrenchment,  channel  substrate  composition, 
and  slope  are  the  stream  reach  morphology 
variables  that  have  the  greatest  influence  on  the 
macroinvertebrate  communities  included  in  this 
study.  Entrenchment  ratio  is  a critical  component 
of  the  Rosgen  stream  classification  system  and  is 
largely  responsible  for  the  system’s  effectiveness 
at  partitioning  biological  variation  among 
HQ-CWFs.  Entrenchment  ratio  describes  the 
vertical  containment  of  the  stream  channel  and 
provides  insight  into  how  stream  energy  is 
dissipated  during  high  flows.  For  example,  the 
top-of-bank  elevations  of  entrenched  streams 
(stream  types  F and  G)  are  much  higher  than 
bankflill  stage  (Figure  13).  At  flows  greater  than 
bankfull  stage,  streamflow  is  contained  within  the 
streambanks  and  stream  depth  increases  much 
faster  than  width.  During  floods,  this  containment 
of  streamflow,  which  results  in  stream  energy 


being  exerted  on  both  streambanks,  is  partly 
responsible  for  the  high  bank  erosion  rates 
observed  in  entrenched  streams  (Rosgen,  1996). 
In  contrast  to  entrenched  streams,  slightly 
entrenched  streams  (stream  types  C and  E) 
overtop  their  streambanks  at  flows  greater  than 
bankfull  stage,  and  stream  width  increases 
substantially  and  stream  energy  is  dissipated  as 
streamflow  moves  onto  the  stream’s  relatively 
broad  floodplain.  The  vertical  containment  of 
moderately  entrenched  streams  (type  B streams)  is 
intermediate  to  that  of  entrenched  and  slightly 
entrenched  streams. 

Entrenchment  ratio,  which  is  the  primary 
stream  reach  morphology  variable  that 
distinguishes  Rosgen  level  I type  B streams  from 
type  C streams,  is  significantly  correlated  to 
taxonomic  richness  and  modified  EPT  index 
scores,  and  taxonomic  richness  is  significantly 
different  between  level  I type  B and  C streams. 
Although  there  is  no  statistically  significant 
difference  between  the  mean  modified  EPT  index 
scores  of  level  I or  modified  level  n stream  types, 
modified  level  n B3  streams  tend  to  have  higher 
modified  EPT  index  scores  than  the  other  stream 
types  (Figure  1 1).  Furthermore,  the  influence  of 
entrenchment  on  EPT  communities  is  reflected  in 
the  significantly  different  modified  percent 
mayfly  metric  scores  of  B and  C streams  (Figure 
9).  The  most  noticeable  difference  in  the  mayfly 
communities  of  B vs.  C streams  is  the  occurrence 
of  the  pollution  intolerant  genus  Paraleptophlebia 
(Ephemeroptera:  Leptophlebiidae).  This  genus  is 
common  in  B streams  (present  at  87  percent  of  the 
B streams),  but  is  relatively  rare  in  type  C streams 
(present  at  23  percent  of  the  C streams). 

The  Rosgen  classification  system  is  based  on 
a suite  of  parameters  that  describe  stream  reach 
channel  morphology  and  floodplain 
characteristics.  For  each  of  these  parameters,  a 
range  of  values  (category  intervals)  were 
developed  based  on  frequency  histograms  of  data 
collected  at  stable,  natural  stream  reaches 
(Rosgen,  1996).  These  category  intervals,  which 
are  broad  enough  to  account  for  the  natural 
variability  in  morphological  conditions  associated 
with  channel  adjustment  processes,  are  used  in  the 
Rosgen  classification  process  to  group  streams 
into  consistent,  reproducible  natural  stream  type 
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Figure  13.  Examples  of  Entrenched,  Moderately  Entrenched,  and  Slightly  Entrenched  Stream  Types 


categories.  Therefore,  in  addition  to  its  utility  as  a 
method  of  identifying  similar  stream  reaches, 
Rosgen  classification  information  is  useful  for 
identifying  streams  that  have  deviated  from  stable, 
natural  morphological  conditions. 

The  Rosgen  classification  data  collected  in 
this  study  indicate  that  the  morphological 
conditions  of  several  sample  sites  are  in 
disequilibrium  and  have  deviated  from  stable, 
natural  conditions.  These  sites  include  Letort 
Spring  Run  (LTRT),  both  sites  on  Hammer  Creek 
(HAMU  and  HAML),  Horse  Hollow  Run 
(HORS),  and  Unnamed  Tributary  to  East  Branch 
Codorus  Creek  (UEBC).  For  example,  LTRT  and 
HAMU  do  not  fit  into  any  of  the  Rosgen 
classification  categories,  but  these  sites  closely 
resemble  G stream  types,  which  Rosgen  (1996) 
describes  as  occurring  in  valleys  such  as  those  of 
LTRT  and  HAMU  under  conditions  of  “instability 
or  disequilibrium  that  are  often  imposed  by 
watershed  changes  and/or  direct  channel  impacts.” 
Several  miles  downstream  of  HAMU,  HAML  also 
is  in  disequilibrium.  Although  HAML  fits  into 
the  F4  category,  it  is  located  in  a type  of  valley 
that  is  typically  drained  by  type  B streams  rather 
than  type  F streams,  indicating  that  this  stream  is 
not  in  equilibrium  with  its  valley.  F stream  types 
occur  naturally  in  canyons,  gorges,  fault- 
controlled  valleys,  glacial  outwash,  and  “U” 
shaped  glacial  troughs. 

With  the  exception  of  Willow  Creek  (WILL) 
and  Hutton  Run  (HUTT),  most  of  the  outliers  in 
the  original  data  set  of  35  sample  sites  show  some 
degree  of  departure  from  stable,  natural 
morphological  conditions.  Willow  Creek,  which 
was  determined  to  be  an  outlier  based  solely  on 
limestone-influenced  water  chemistry,  is  the  only 

limestone  influenced  site  with  stable,  natural 

morphological  characteristics.  In  addition  to 
stable  morphological  conditions,  the 

macroinvertebrate  community  of  Willow  Creek  is 
characterized  by  metric  scores  comparable  to 
other  type  C streams,  with  the  exception  of  a very 
high  modified  percent  mayfly  index  score  due  to  a 
high  abundance  of  the  genus  Isonychia 
(Ephemeroptera;  Isonychiidae).  This  relationship 
observed  between  stream  reach  morphological 
conditions  and  the  biological  integrity  of  the 
limestone  influenced  streams  included  in  this 


study  (WILL  vs.  LTRT,  HAMU,  and  HAML) 
suggests  that  Rosgen  classification  data  contains 
valuable  information  that  could  enhance 
evaluations  of  the  ecological  significance  of 
streams  and  rivers. 

The  Rosgen  level  I delineations  referred  to  in 
this  study  were  determined  based  on  field  data 
gathered  during  level  II  assessments.  However, 
level  I classification  is  designed  to  be  a coarse- 
scale  geomorphic  characterization  that  can  be 
performed  using  information  taken  from 
topographic  and  landform  maps  and  aerial 
photographs,  followed  by  subsequent  verification 
with  field  data  from  selected  reference  reaches. 
Thus,  level  I classifications  performed  in  the 
office  and  verified  in  the  field  may  provide  water 
resource  managers  with  an  accurate,  cost  effective 
procedure  for  classifying  streams  and  screening 
stream  reach  morphological  conditions. 
Furthermore,  the  efficiency  of  Rosgen  level  I 
classification  and  data  management  could  be 
greatly  enhanced  if  these  activities  can  be 
successfully  conducted  using  GIS  and  spatial 
analysis  software. 

MANAGEMENT  IMPLICATIONS 

The  data  included  in  this  study  were  useful  for 
identifying  several  relationships  between  stream 
reach  morphology  and  the  macroinvertebrate 
communities  of  HQ-CWFs.  Channel  substrate 
characteristics,  entrenchment,  and  surface  water 
slope  are  the  stream  reach  morphology  features 
observed  to  have  the  greatest  influence  on  the 
metric  scores  of  these  streams.  Taxa  richness, 
modified  EPT,  and  modified  Hilsenhoff  index 
scores  are  influenced  by  channel  substrate 
characteristics.  Taxa  richness  and  modified  EPT 
index  scores  are  influenced  by  entrenchment  and 
modified  Hilsenhoff  index  scores  are  influenced 
by  surface  water  slope. 

The  Rosgen  Stream  Classification  System 
provides  a means  for  classifying  and  assessing 
streams  based  on  stream  reach  morphology 
variables  that  substantially  influence  the  physical 
characteristics  and  benthic  macroinvertebrate 
communities  of  HQ-CWFs.  Whereas  no 
significant  relationships  exist  between  stream 
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order  and  metric  scores,  Rosgen  classification 
designation  is  superior  to  stream  order  designation 
as  a method  for  partitioning  biological  variability 
in  HQ-CWFs  and  classifying  these  streams  at  the 
reach  spatial  scale. 

The  relationships  observed  between  stream 
reach  morphology  and  the  biological  integrity’  of 
the  streams  included  in  this  study,  limestone 
influenced  streams  in  particular,  indicate  that 
Rosgen  classification  data  contains  valuable 
information  pertaining  to  the  general, 
morphological  “state”  or  condition  of  stream 
reaches.  This  information  provides  objective, 
quantitative  data  for  documenting  the 
morphological  conditions  of  streams  with  regard 
to  deviations  fi-om  stable,  natural  morphological 
conditions.  The  potential  inclusion  of  this 
information  in  bioassessments  and/or  ecological 
significance  evaluations  deserves  further 
investigation. 

Other  potential  uses  of  Rosgen  classification 
data  include  applications  related  to  hydraulic 
geometry,  flow  resistance,  grazing  management 
techniques,  fish  habitat  improvement  structure 
suitability,  stream  channel  restoration  design,  and 
watershed  management  activities  based  on  stream 
type.  Of  particular  interest,  is  the  potential  utility’ 
of  Rosgen  stream  classification  data  to  develop 
watershed  management  guidelines.  Myers  and 
Swanson  (1992)  identified  how  different  Rosgen 
stream  types  responded  to  riparian  grazing. 
Rosgen  (1994)  expanded  on  the  findings  of 
Meyers  and  Swanson  and  summarized  the  affects 
of  grazing  on  different  stream  types  with  respect 
to  their  sensitivity  to  disturbance,  recovery 
potential,  sediment  supply,  streambank  erosion 
potential,  and  the  influence  of  vegetation  at 
controlling  stability. 

In  categorizing  the  riparian  grazing  data  by 
stream  type,  Rosgen  developed  a framework  for 
identifying  stream  reaches  that  are  most  sensitive 
to  disturbance  and  have  the  highest  recovery 
potential,  based  on  stream  type.  It  is  highly 
probable  that,  if  the  relationships  between 
additional  watershed  activities  (e.g.,  land  use/land 
cover  changes,  riparian  vegetation  conditions,  the 
occurrence  of  specific  farming  practices)  and 
Rosgen  stream  type  designation  were  identified, 
stream  classification  could  be  used  as  a 


framework  for  organizmg  watershed  data, 
developing  watershed  management  guidelines, 
and  effectively  targeting  a variety'  of  watershed 
management  and  restoration  activities. 

Future  investigation  comparing  the  Rosgen 
classification  data  obtained  in  this  study  w'ith 
classification  data  for  the  same  stream  reaches 
generated  in  the  office  would  be  valuable  for 
determining  whether  office  delineations,  possibly 
using  GIS  software,  are  accurate  for  classitying 
stream  reaches  and  describing  stream  reach 
morphology.  If  so,  GIS-generated  Rosgen 
classification  information  could  provide  water 
resource  managers  with  a cost-effective  procedure 
for  classifying  streams,  screening  stream  reach 
morphological  conditions,  organizing  watershed 
data,  and  targeting  watershed  management 
activities. 
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Appendix  A 

Channel  Cross-Section  and  Pebble  Count  Data 
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Figure  Al.  Cross-Section  and  Pebble  Count  Data  From  Asaph  Run 
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Figure  A2.  Cross-Section  and  Pebble  Count  Data  From  Baldwin  Branch 
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Figure  A3,  Cross-Section  and  Pebble  Count  Data  From  Cork  Hill  Creek 
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Figure  A4.  Cross-Section  and  Pebble  Count  Data  From  Deep  Hollow 
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Figure  A5.  Cross-Section  and  Pebble  Count  Data  From  Dead  Woman  Flollow 
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Figure  A6.  Cross-Section  and  Pebble  Count  Data  From  East  Branch  Codorus  Creek 
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Figure  A 7.  Cross-Section  and  Pebble  Count  Data  From  East  Branch  Lackawanna  River 
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Figure  A8.  Cross-Section  and  Pebble  Count  Data  From  East  Branch  Roaring  Brook 
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Figure  A9.  Cross-Section  and  Pebble  Count  Data  From  Gazzam  Run 
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Figure  AlO.  Cross-Section  and  Pebble  Count  Data  From  Greenlee  Run 
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Figure  AIL  Cross-Section  and  Pebble  Count  Data  From  Green  Run  (Lower) 
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Figure  A 12.  Cross-Section  and  Pebble  Count  Data  From  Green  Run  (Upper) 
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Figure  A1 3.  Cross-Section  and  Pebble  Count  Data  From  Hammer  Creek  (Lower) 
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Figure  A14.  Cross-Section  and  Pebble  Count  Data  From  Hammer  Creek  (Upper) 
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Figure  A1 5.  Cross-Section  and  Pebble  Count  Data  From  Harveys  Creek 


51 


NUMBER  OF  PARTICLES 


CUMULATIVE  PERCENT  FINER  ELEVATION  (ft) 


CROSS-SECTION  DATA 


PEBBLE-COUNT  DATA 


Figure  A 16.  Cross-Section  and  Pebble  Count  Data  From  Horse  Hollow 
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Figure  A 17.  Cross-Section  and  Pebble  Count  Data  From  Hutton  Run 
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Figure  A1 8.  Cross-Section  and  Pebble  Count  Data  From  Iron  Run 
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Figure  A 19.  Cross-Section  and  Pebble  Count  Data  From  Lebo  Run 
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Figure  A20.  Cross-Section  and  Pebble  Count  Data  From  Little  Moores  Run 
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Figure  All.  Cross-Section  and  Pebble  Count  Data  From  Letort  Spring  Run 
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Figure  A22.  Cross-Section  and  Pebble  Count  Data  From  Mountain  Creek 
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Figure  A23.  Cross-Section  and  Pebble  Count  Data  From  Murray  Run 
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Figure  A24.  Cross-Section  and  Pebble  Count  Data  From  Pequea  Creek 
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Figure  A25.  Cross-Section  and  Pebble  Count  Data  From  Rice  Run 
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Figure  A26.  Cross-Section  and  Pebble  Count  Data  From  Sixmile  Run 
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Figure  A27,  Cross-Section  and  Pebble  Count  Data  From  Standing  Stone  Creek 
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Figure  A28.  Cross-Section  and  Pebble  Count  Data  From  Tipton  Run 
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Figure  A29.  Cross-Section  and  Pebble  Count  Data  From  Toms  Run  (Lower) 
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Figure  A30.  Cross-Section  and  Pebble  Count  Data  From  Toms  Run  (Upper) 
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Figure  A3L  Cross-Section  and  Pebble  Count  Data  From  Unnamed  Tributary  to  East  Branch 
Codorus  Creek 
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Figure  A32.  Cross-Section  and  Pebble  Count  Data  From  Unnamed  Tributary  to  Harveys  Creek 
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Figure  A33.  Cross-Section  and  Pebble  Count  Data  From  Unnamed  Tributary  to  Roaring  Brook 
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Figure  A34.  Cross-Section  and  Pebble  Count  Data  From  Willow  Run 
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Figure  A35.  Cross-Section  and  Pebble  Count  Data  From  White  Deer  Creek 
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Leucrocuta 

Rhithrogena 

Stenacron 

Stenonema 

6 

15 

15 

3 

8 

Isonychiidae 

Isonychia 

8 

6 

Leptophlebiidae 

Paraleptophlebia 

8 

5 

6 

Tricorythidae 

Tricorythodes 

Megaloptera 

Corydalidae 

Nigronia 

3 

1 

6 

2 

2 

SiaUdae 

Sialis 

1 

Odonata 

Aeshnidae 

Boyeha 

Gomphidae 

Lanthus 

1 

2 

1 

Stylogomphus  albistylus 

1 

Plecoptera 

Chloroperlidae 

Alloperla 

9 

Haploperla 

Sweltsa 

Leuctridae 

Leuctra 

Peltoperiidae 

Peltoperla 

Tallaperla 

4 

Perlidae 

Acroneuria 

7 

1 

6 

3 

8 

Agnetina 

5 

Eccoptura 

4 

1 

Paragnetina 

Perlodidae 

Cultus 

Diploperla 

Hydroperla 

Isoperla 

1 

1 

2 

1 

Remenus 

2 

Pteronarcyidae 

Pteronarcys 

1 

1 

T aeniopterygidae 

Sv-ophopteryx 

Taemopteryx 

6 

10 

1 
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Class:  Order 

Family 

Family/Genus 

DWMN 

EBCC 

EBLR 

EBRB 

GAZZ 

GRLE 

Trichoptera 

Brachycentridae 

Brachycentrus 

Glossosomatidae 

As,apetus 

1 

11 

Glossosoma 

1 

Protoptila 

3 

Hydropsychidae 

Cheumatopsyche 

4 

13 

2 

1 

7 

2 

Diplectrona 

3 

8 

2 

HomoplecP-a 

Hydropsyche 

30 

21 

8 

9 

22 

Potamyia  flava 

2 

25 

Lepidostomatidae 

Theliopsyche 

Odontoceridae 

Psilop-eta 

Philopotamidae 

Chimarra 

1 

Dolophilodes 

2 

2 

7 

6 

4 

25 

Polycentropodidae 

Neureclipsis 

3 

1 

PolycenPopus 

Psychomyiidae 

Lype  diverse 

Rhyacophilidae 

Rhyacophila 

1 

Rliynchobdellida 

Glossiphonidae 

BaPachobdella 

Crustacea:  Amphipoda 

Gamma  ridae 

Gammarus 

Decapoda 

Cambaridae 

Cambarus 

2 

Isopoda 

Asellidae 

Caecidotea 

Lirceus 

Gastropoda:  Gastropoda 

Pleuroceridae 

Leptoxis 
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Class:  Order  | 

Family  | 

Family/Genus 

GRNL 

GRNU 

HAML 

HAMU 

HARV 

HORS 

Insecta:  Coleoptera 

Clmidae 

Optioservus 

1 

3 

22 

8 

6 

Oulimnius 

°romoresia 

1 

Stenelmis 

4 

3 

Psephenidae 

Psephenus 

Ectopria 

Ptilodactylidae 

Anchytarsus 

Diptera 

Athericidae 

Athenx 

Ceratopoeonidae 

Bezzia 

Chironomidae 

Chironomidae 

7 

6 

3 

34 

11 

Empididae 

Hemerodromia 

Simuliidae 

Simuliidae 

1 

2 

2 

Tabanidae 

Tabanus 

1 

Tipulidae 

Antocha 

4 

Dicranota 

Hexatoma 

1 

1 

2 

Limonia 

Tipula 

1 

Ephemeroptera 

Ameletidae 

Ameletus 

1 

Baetidae 

Acentrella 

3 

Baetis 

12 

8 

15 

Centroptilum 

Cloeon 

Baetiscidae 

Baetisca 

Ephemerellidae 

Ephemerella 

Eurylophella 

2 

Serratella 

3 

1 

Ephemeridae 

Ephemera 

1 

Heptaeeniidae 

Epeorus 

Heptazenia 

Lzucrocuta 

Rhithrozena 

Slenacron 

Stenonema 

4 

13 

10 

21 

Isonychiidae 

hortychia 

17 

Leptophlebiidae 

Paraleptophleb  ia 

8 

6 

1 

Tricorythidae 

Tricorythodes 

1 

Megaloptera 

Corydalidae 

Nizronia 

1 

2 

1 

9 

Sialidae 

Sialis 

Odonata 

Aeshnidae 

Boyeria 

1 

2 

C^mphidae 

Lanthus 

1 

Stylozomphus  albistylus 

1 

Plecoptera 

Chloroperlidae 

AUoperla 

2 

Haploperla 

Sweltsa 

Lenctridae 

LeucP-a 

Peltoperildae 

Peltoperla 

TaUaperla 

1 

Perlidae 

Acroneuria 

6 

Aznetina 

Eccoptura 

2 

Paraznetina 

10 

Perlodidae 

Cultus 

Diploperla 

3 

Hydroperla 

Isoperla 

2 

Remenus 

Pteronarcyidae 

Pteronarcys 

T aeniopterygidae 

SP'ophopteryx 

Taeniopteryx 

3 
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Class:  Order 

Family 

Family/Genus 

GRNL 

GRNU 

HAML 

HAMU 

HARV 

HORS 

Trichoptera 

Brachycentridae 

Brachycentrus 

Glossosomatidae 

Agapetus 

Glossosoma 

Protoptila 

Hydropsychidae 

Cheumatopsyche 

17 

15 

9 

Diplectrona 

4 

4 

1 

Homoplectra 

Hydropsyche 

6 

3 

45 

23 

6 

Potamyia  flava 

8 

16 

Lepidostomaddae 

Theliopsyche 

Odontoceridae 

Psilotreta 

Philopotamidae 

Chimarra 

1 

21 

Dolophilodes 

15 

11 

Poiycentropodidae 

Neureclipsts 

Polycentropus 

Psychomyildae 

Lype  di versa 

2 

Rhyacophilidae 

Rhyacophila 

1 

4 

Rhynchobdellida 

Glossiphonidae 

Batrachobdella 

Crustacea:  Amphipoda 

Gammaridae 

Gammarus 

18 

1 

Decapoda 

Cambaridae 

Cambarus 

1 

2 

Isopoda 

Asellidae 

Caecidotea 

2 

Lirceus 

Gastropoda:  Gastropoda 

Pleuroceridae 

Leptoxis 

86 


Class:  Order 

Family 

Family/Genus 

HUTT 

IRON  1 LEBO 

LMOR 

LTRT 

MNTN 

Insecta:  Coleoptera 

ilinidae 

Optioservus 

2 

1 

6 

Dulimmus 

3 

1 

Promoresia 

4 

Stenelmis 

1 

Psephenidae 

Psephenus 

Ectopria 

2 

1 

Ptilodactylidae 

Anchytarsus 

Diptera 

Athericidae 

Atherix 

20 

1 

Ceratopogonidae 

Be::zia 

Chironcmldae 

Chironomidae 

1 

3 

2 

7 

2 

15 

Empididae 

Hemerodromia 

Simuliidae 

Simuliidae 

1 

1 

Tabanidae 

Tabanus 

Tipulidae 

Antocha 

Dicranota 

2 

1 

2 

Hexatoma 

9 

5 

10 

Limonia 

Tipula 

3 

1 

1 

1 

1 

Ephemeroptera 

Ameletidae 

Ameletus 

Baetidae 

Acentrella 

Baetis 

26 

8 

2 

1 

Centroptilum 

Cloeon 

Baetiscidae 

Baetisca 

Ephemerellidae 

Ephemerella 

15 

1 

1 

Eurylophella 

Serratella 

4 

Ephemeridae 

Ephemera 

Heptageniidae 

Epeorus 

3 

Hepta^enia 

3 

13 

Leucrocuta 

4 

Rhithrogena 

Stenacron 

Stenonema 

3 

1 

4 

8 

Isonychiidae 

honychia 

L«ptophlebiidae 

Paraleptophlebia 

1 

T ricorythidae 

Tricorythodes 

Megaloptera 

Corydalidae 

Nigronia 

2 

1 

Sialidae 

Sialis 

Odonata 

Aeshnidae 

Boyeria 

Gomphidae 

Lanthus 

5 

Stylogomphus  albistylus 

Plecoptera 

ChloroperUdae 

Alloperla 

Haploperla 

Sweltsa 

Leuctridae 

Leuctra 

Peltoperlidae 

Peltoperla 

Tallaperla 

3 

3 

8 

Perlidae 

Acroneurw 

2 

Agnetina 

2 

Eccoptura 

6 

1 

Paragnetina 

Periodldae 

Cultus 

Diploperla 

Hydroperla 

2 

4 

1 

Isoperla 

6 

1 

3 

4 

Remenus 

Pteronarcyidae 

Pteronarcys 

3 

4 

1 

T aenioptery  gidae 

Stroohopteryx 

1 

Taemopteryx 

87 


Class:  Order 

Family 

Family/Genus 

HUTT 

IRON 

LEBO 

LMOR 

LTRT 

MNTN 

Trichoptera 

Brachycentridae 

Brachvcentrus 

1 

1 

Glossosomatidae 

A^apetus 

Glossosoma 

2 

1 

Protoptila 

Hydropsychidae 

Cheumatopsyche 

17 

Diplectrona 

11 

1 

3 

Homoplecp-a 

Hydropsyche 

15 

4 

9 

10 

Potamyia  flava 

15 

6 

S 

Lepidostomatidae 

Theliopsyche 

6 

Odontoceridae 

Psilotreta 

Philopotamidae 

Chimarra 

Dolophilodes 

10 

13 

5 

2 

Polycentropodidae 

Neureclipsis 

Polycentropus 

Psychomyiidae 

Lype  diversa 

Rhyacophilidae 

Rhyacophila 

2 

3 

1 

Rhynchobdellida 

Glossiphonidae 

Batrachobdella 

Crustacea:  Amphipoda 

Gammaiidae 

Gammarus 

63 

Decapoda 

Cambaridae 

Cambarus 

1 

1 

1 

Isopoda 

Asellidae 

Caecidotea 

Lirceus 

37 

Gastropoda:  Gastropoda 

Pleuroceridae 

Leptoxis 
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Class:  Order 

Familv 

Family/Genus  | MURY 

PEQU 

RICE 

SIXM 

STSN 

TIPT 

Insecta:  Coleoptera 

iJmidae 

Dptioservus 

1 1 3 

( 

Oulimnius 

1 

1 

^romoresia 

2 

Stenelmis 

^sephenidae 

Psevhenus 

1 

1 

j 

Ectopna 

2 

3 

Ptilodactylidae 

Anchytarsus 

Diptera 

Atheriddae 

Atherix 

Ceratopogonidae 

Bezzia 

Chironomidae 

Chironomidae 

12 

1 

3 

5 

6 

8 

Empididae 

Hemerodromia 

1 

Simuliidae 

Simuliidae 

13 

4 

1 

Tabanidae 

Tabanus 

Tipnlidae 

Antocha 

6 

Dicranota 

Hexatoma 

2 

6 

2 

1 

Limonia 

Tipula 

2 

Ephemeroptera 

Ameletidae 

Ameletus 

Baetidae 

Acennella 

1 

2 

Baetis 

4 

2 

21 

CentToptilum 

Cloeon 

4 

Baetisddae 

Baetisca 

1 

Ephemerellidae 

Ephemerella 

2 

Eurylophella 

Serratella 

1 

Ephemeridae 

Ephemera 

Heptageniidae 

Epeorus 

Hepta^enia 

Leucrocuta 

4 

Rhithrogena 

Stenacron 

Stenonema 

13 

4 

8 

2 

3 

Isonychiidae 

Isortychia 

1 

3 

L«ptophlebiidae 

Paraleptophlebia 

2 

1 

1 

Tricorythidae 

Tricorythodes 

Megaloptera 

Corydalidae 

Nigronia 

1 

5 

3 

1 

Sialidae 

Sialis 

1 

Odonata 

Aeshnidae 

Boyeria 

Gomphidae 

Lanthus 

4 

1 

Stylogomphus  albistylus 

1 

1 

Plecoptera 

ChloroperUdae 

Alloperla 

4 

4 

Haploperla 

Sweltsa 

Leuctridae 

Leucp-a 

4 

Peltoperlidae 

Peltoperla 

2 

Tallaperla 

1 

1 

2 

Periidae 

Acroneuria 

7 

8 

Agnetjna 

Eccoptura 

9 

1 

Paragnetina 

2 

Perlodldae 

Cultus 

Diploperla 

Hydroperla 

Isoperla 

9 

5 

Remenus 

2 

Pteronarcyidae 

Pteronarcys 

19 

3 

1 

T aenioptery  gidae 

Strophopteryx 

Taemopteryx 
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Class:  Order 

Family 

Family/Genus 

MURY 

PEQU 

RICE 

SIXM 

STSN 

TIPT 

Trichoptera 

Brachycentridae 

Brachycentrus 

Glossosomatidae 

A^apetus 

Glossosoma 

1 

Protoptila 

Hydropsychidae 

Cheumatopsyche 

1 

2 

4 

1 

14 

DiplecP-ona 

3 

1 

9 

Homoplectra 

Hydropsyche 

46 

19 

4 

6 

7 

Potamyia  flava 

3 

7 

3 

12 

3 

Lepidostomatidae 

Theliopsyche 

Odontoceridae 

Psilotreta 

Philopotamidae 

Chimarra 

6 

Dolophilodes 

21 

4 

5 

21 

24 

41 

Polycentropodidae 

Neureclipsis 

Polycentropus 

2 

Psychomyiidae 

Lype  diverse 

2 

Rhyacophilidae 

Rhyacophila 

3 

4 

10 

10 

3 

Rhynchobdellida 

Glossiphonidae 

Batrachobdella 

1 

Crustacea;  Amphipoda 

Gammaridae 

Gammarus 

Decapoda 

Cambaridae 

Cambarus 

1 

1 

Isopoda 

Asellidae 

Caecidotea 

Lirceus 

Gastropoda:  Gastropoda 

Pleuroceridae 

Leptoxis 
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Class:  Order  | Family 

Family/Genus 

TOML 

TOMU 

UEBC 

UHAR 

URRB 

WILL 

WTDR 

Insecta:  Coleoptera 

Clmidae 

Optioservus 

4 

17 

1 

10 

Oulimnius 

2 

3 

Promoresia 

4 

1 

10 

Stenelrms 

Psephenidae 

Psephenus 

1 

Ectopha 

2 

Ptilodactylidae 

Anchytarsus 

1 

1 

1 

Diptera 

Athericidae 

Atherix 

10 

Ceratopogonidae 

Bezzia 

Chironomidae 

Chironomidae 

4 

4 

15 

6 

9 

Empididae 

Hemerodromia 

Simuliidae 

Simuliidae 

3 

18 

Tabanidae 

Tabanus 

Tipalidae 

Antocha 

1 

Dicranota 

6 

1 

1 

8 

Hexatoma 

2 

3 

3 

3 

Limonia 

1 

Tipula 

6 

1 

Ephemeroptera 

Ameletidae 

Ameletus 

Baetidae 

Acenp-ella 

Baetis 

2 

2 

Centroptilum 

2 

Cloeon 

2 

Baetiscidae 

Baetisca 

Ephemerellidae 

Ephemerella 

7 

4 

Eurylophella 

Serratella 

Ephemeridae 

Ephemera 

1 

Heptageniidae 

Epeorus 

1 

1 

Heptagenia 

Leucrocuta 

Rhithro^ena 

Stenacron 

1 

Stenonema 

1 

5 

19 

Isonychiidae 

Isonychia 

30 

Leptophlebiidae 

Paraleptophlebia 

21 

5 

2 

Tricorythidae 

Tricorythodes 

Megaloptera 

Corydalidae 

Nigronia 

2 

6 

2 

1 

5 

9 

Sialidae 

Sialis 

Odonata 

Aeshnidae 

Boyeria 

1 

Gomphidae 

Lanthus 

2 

Stylo^omphus  albistyhis 

1 

Plecoptera 

Chloroperiidae 

Alloperla 

2 

Haploperla 

Sweltsa 

Leuctridae 

Leuctra 

2 

PeWoperiidae 

Peltoperla 

Tallaperla 

1 

3 

Perlldae 

Acroneuria 

1 

3 

6 

12 

1 

3 

Aznetina 

Eccoptura 

5 

3 

Para^netina 

Perlodidae 

Cultus 

1 

Diploperla 

Hydroperla 

6 

Isoperla 

1 

Remenus 

Pteronarcyidae 

Pteronarcys 

1 

1 

2 

T aenioptery  gidae 

SP-ophopteryx 

1 

7 

Taenwpteryx 
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Class:  Order 

Family 

Family/Genus 

TOML 

TOMU 

UEBC 

UHAR  1 URRB  WILL 

WTDR 

T richoptera 

Brachycentridae 

Brachycentrus 

4 

Glossosomatidae 

A^apetus 

Glossosoma 

4 

5 

1 

Protoptila 

Hydropsychidae 

Cheumatopsyche 

4 

4 

7 

14 

7 

5 

9 

Diplectrona 

3 

1 

2 

2 

Homoplectra 

4 

Hydropsyche 

2 

23 

5 

3 

3 

4 

Potamyia  flava 

4 

11 

2 

Lepidostomatidae 

Theliopsyche 

1 

Odontoceridae 

Psilotreta 

3 

Philopotaniidae 

Chimarra 

19 

1 

4 

Dolophilodes 

2 

1 

13 

11 

16 

2 

16 

Polycentropodidae 

Neureclipsis 

Polycentropus 

Psychomyiidae 

Lype  diverse 

3 

Rhyacophilidae 

Rhyacophila 

1 

2 

1 

2 

5 

Rhynchobdellida 

Glossiphonidae 

Batrachobdella 

1 

Crustacea:  Amphipoda 

Ganunaridae 

Gammarus 

Decapoda 

Cambaridae 

Cambarus 

2 

1 

Isopoda 

Asellidae 

Caecidotea 

Lirceus 

Gastropoda:  Gastropoda 

Pleuroceridae 

Leptoxis 

1 
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